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^r\' Abstract 

In a recent paper [Tj, we have studied the vacuum solutions of 
Bianchi types / and V spacetimes in the framework of metric f(R) 
qq ■ gravity. Here we extend this work to perfect fluid solutions. For this 

0^ . purpose, we take stiff matter to find energy density and pressure of the 



universe. In particular, we find two exact solutions in each case which 

correspond to two models of the universe. The first solution gives 
f ^ ■ a singular model while the second solution provides a non-singular 

^^ \ model. The physical behavior of these models has been discussed 

using some physical quantities. Also, the function of the Ricci scalar 

is evaluated. 

X 
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1 Introduction 

The most interesting phenomenon, which the modern day physics deals with 
today, is the accelerated expansion of the universe. The explanation for this 
universe expansion has both theoretical as well as experimental background 



*msharif@math.pu.edu.pk 
tfrasat@hotmail.com 



[2J. It has been found that most of the universe contains dark energy and 
dark matter. The cosmological constant is shown [3] to be an alternative 
for dark energy. f{R) theory of gravity is an attractive candidate as an 
alternative theory of gravity in which a general function of the Ricci scalar, 
f(R), replaces R in the standard Einstein-Hilbert Lagrangian. This provides 
[I] an easy unification of early time inflation and late time acceleration. The 
theory also gives a natural gravitational alternative to dark energy. 

Lobo [5] gave a brief review of some of the modified theories of gravity 
that address dark energy and the dark matter problems. In another review 
[6], Faraoni discussed the explanation of the cosmic acceleration alternative 
to dark energy in the various versions of f(R) theories of gravity. He inves- 
tigated the successes of f(R) gravity together with the challenges imposed 
by minimal criteria for their viability. Sotiriou and Faraoni [7j presented 
some important aspects of f(R) theories of gravity in Metric, Palatini and 
Metric- Affine formalisms. They discussed the motivation, actions, field equa- 
tions, equivalence with other theories, cosmological aspects and constraints, 
viability criteria and astrophysical applications. 

The f(R) theory of gravity is one of the modified theories which is con- 
sidered most suitable due to cosmologically important f(R) models. These 
models consist of higher order curvature invariants as functions of the Ricci 
scalar. Viable f(R) gravity models [8] have been proposed which show the 
unification of early-time inflation and late-time acceleration. It is hoped that 
the problem of dark matter can be addressed by using viable f(R) gravity 
models. Multamaki and Vilja [9] investigated spherically symmetric vacuum 
solutions in f(R) theory. The same authors [10] also studied the perfect 
fluid solutions and showed that pressure and density did not uniquely deter- 
mine f(R). Cognola et al. [TT] investigated f(R) gravity at one-loop level 
in de-Sitter universe. Capozziello et al. [12] explored spherically symmetric 
solutions of f(R) theories of gravity via the Noether symmetry approach. 
Recently, Sharif and Kausar [13] studied non- vacuum static spherically sym- 
metric solutions in f(R) gravity. 

Hollenstein and Lobo [T3j analyzed exact solutions of static spherically 
symmetric spacetimes in f(R) gravity coupled to non-linear electrodynam- 
ics. Azadi et al. [15] studied cylindrically symmetric vacuum solutions in 
this theory. Momeni and Gholizade [16] extended cylindrically symmetric 
solutions in a more general way. We have explored static plane symmetric 
vacuum solutions [T7] in f(R) gravity. The field equations are solved using 
the assumption of constant scalar curvature which may be zero or non- zero. 



However, very few attempts have been made so far for solutions with non- 
constant scalar curvature. 

Bianchi types I and V spacetimes play an important role in the study 
of universe. Lorenz-Petzold [18] studied exact Brans-Dicke Bianchi type / 
solutions with a cosmo logical constant. The same author [19] also derived 
solutions for locally rotationally symmetric (LRS) Bianchi type V space- 
time in the Brans-Dicke theory of gravitation. It was found that solutions 
represented anisotropic cosmological models filled with stiff matter and an 
electromagnetic null field. Berman [20] introduced a different method to 
solve the field equations by using the variation law of Hubble parameter. 
The main feature of the variation law is that it gives a constant value of 
the deceleration parameter. Kumar and Singh [2[T] investigated perfect fluid 
solutions using Bianchi type / spacetime in scalar-tensor theory. Singh [22] 
investigated LRS Bianchi type V cosmology with heat flow in Scalar- Tensor 
theory. Paul et al. [23] investigated FRW cosmologies in f(R) gravity. In 
a recent paper [1] , we have studied the vacuum solutions of Bianchi types / 
and V spacetimes in the framework of f{R) gravity. 

In this paper, we focuss our attention to explore the perfect fluid solutions 
of Bianchi types / and V spacetimes in metric f{R) gravity. For this purpose, 
we take the case of stiff matter to find energy density and pressure of the 
universe. The paper is organized as follows: In section 2, we briefly give the 
field equations in metric f(R) gravity. Sections 3 and 4 are used to find some 
exact solutions of Bianchi types / and V spacetimes respectively. In the last 
section, we discuss the results. 



2 f(R) Gravity Formalism 

The action for f(R) gravity is given by 



S = J ^9(j^Qf( R ) + L m)d 4 x, (1) 

where f{R) is a general function of the Ricci scalar and L m is the matter 
Lagrangian. The corresponding field equations are the following: 

F(R)R^ - i/(ii)v, - V^V,F( J R) + g^DF(R) = kT^ (2) 

where F(R) = df(R)/dR, □ = V^V^, V^ is the covariant derivative and T^ u 
is the standard matter energy-momentum tensor derived from the Lagrangian 



L m . When we contract the field equations, it follows that 

F(R)R - 2f(R) + 3DF(R) = kT. (3) 

Using this equation in Eq.(j2]), the field equations take the form 

F{R)R, V - V,V U F{R) - kT, v = [ F(R)R - °^ R) - KT }g, u . (4) 

Thus we have eliminated f(R) from the field equations which helps us to 
solve the field equations. 

3 Exact Bianchi Type I Solutions 

In this section we find exact solutions of Bianchi I spacetime in f(R) gravity 
and some physical quantities. 

3.1 Field Equations and Some Physical Quantities 

The line element of Bianchi type I spacetime is given by 

ds 2 = dt 2 - A 2 (t)dx 2 - B 2 (t)dy 2 - C 2 (t)dz 2 , (5) 

where A, B and C are cosmic scale factors. The corresponding Ricci scalar 

is 

n J B C AB BC CA^ 

R =~ 2[ A + B + C + AB + BC + CA^ (6) 

where dot means derivative with respect to t. The energy-momentum tensor 
for perfect fluid yields 

Fy.v = (p + p)u^u v - pg^, (7) 

satisfying the equation of state 

p = wp, < w < 1, (8) 

where p and p are energy density and pressure of the fluid while u^ = 
v /poo(l, 0,0,0) is the four- velocity in co-moving coordinates. Since the met- 
ric ()5]) depends only on t, Eq.(@| yields a set of differential equations for 



F(t), A, B, C, p and p. Thus the subtraction of the 00-component and 
11-component gives 

B C Ab CA AP F k, 
-B-C + AB + CA + AF-F-F {P + P) = °- (9) 

Similarly, we get two more independent equations 

A c Ab bc bf f k, 

-A'C + AB + BC + BF-F-F i ^''> = - (10) 

A B BC CA CF F k, , „ 
-A--B + W + C7A + CF-F-F i " + P) = - (U) 

The conservation equation, T^ u u = 0, leads to 

P+(P + P)[j + B + C\ = °- ( 12 ) 

Finally, we have four differential equations with six unknowns namely A, B, 
C, F, p and p. The solution of these equations is discussed in the next 
subsection. In the following we give definition of some physical quantities. 
We define the average scale factor and the volume scale factor respectively 



a = VABC, V = a 3 = ABC. (13) 

The generalized mean Hubble parameter H is given in the form 

H = ^(H 1 + H 2 + H 3 ), (14) 

where Hi — 4, H 2 = -§, H 3 — Q are the directional Hubble parameters in 
the directions of x, y and z axis respectively. Using Eqs.f ll3|) and (114]) . we 
obtain 

H=~= 1 -(H 1 + H 2 + H 3 ) = -. (15) 

3V 3 a 

The mean anisotropy parameter A is given by 
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^E(^V^) 2 - ( 16 ) 



3^- v H 



The expansion scalar 9 and shear scalar a 2 are defined as follows 

e = < = ^ + l + §- (17) 



1 .... l rr A,o ,B s0 ,C\ AB BC CA 

a" = 

where 



s = r^^-^^B^^-'-XB-w-^ (18) 



<V = T^wK + u v . a h a ^) - -9h^, (19) 

h^v = g^u — u^Uy is the projection tensor. In thermodynamics, the entropy 
of the universe is given by 

TdS = d(pV) + pdV. (20) 

3.2 Solution of the Field Equations 



Subtracting Eqs.pDjI, (01} and (01} from Eqs.©, (05]) and (Ej), we get re- 
spectively 

A B C .A B^ F.A B. n 

B C A.B C\ F.B C\ . , 

B-C + I ( B-C )+ F ( B-C )= °- (22) 

i C B.A C\ F.A c\ . , 

A-C + B ( I-C )+ F ( I-C )=0 - (23) 

These equations are exactly the same as given by Eqs.(18)-(20) in pp. Thus 
we can write the metric functions explicitly as 

f dt 

A = ap 1 exp[q 1 j —}, (24) 

a?F^ ( 25 ) 

f dt 
C = ap 3 exp[q 3 / — ], (26) 

where 

Pi = (dr 2 d2 _1 )5, p 2 = (di^ -1 )^, P3 = (did2 2 )', P\P2P-i = l (27) 
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and 

2ci + c 2 ci - c 2 ci + 2c 2 , . 

<h = g — , 92 = — g — , Q3 = — g — , gi + g2 + g3 = o, (28) 

Q and dj are constants of integration. Using power law relation between F 
and a [I], we have 

F = ka m , (29) 

where k is the constant of proportionality, m is any integer (here taken as 2) 
and a is given by 

a = (nlt + kx)*, n^O 

a = fc 2 exp(/t), n = 0, (30) 

fci and fc 2 are constants of integration. It is mentioned here that we have 
used H = la~ n , I > 0, n > to get the above equation. Thus we obtain 
two values of the average scale factor corresponding to two different models 
of the universe. 

3.3 Model of the Universe when n^O 

For this model, F becomes F = k(nlt + fci)~» and the corresponding metric 
coefficients A, B and C turn out to be 

a ( u , 7 ^ rqi(nlt + fa) 2 ^ 

A = pi{nlt + Ki)« exp[- 



kl(n-l) 

.q 2 (nlt + /ci) 1 ^ - 

kl(n- 1) 

r- i n i / ^ ^3(^ + ^1)^ 

C = Psynlt + A4J" exp[ 



B = P2{nlt + ki) n exp 



7^1 


(31) 


7^1 


(32) 


7^1- 


(33) 



kl(n-l) J 
The mean generalized Hubble parameter and the volume scale factor become 

^SmW "=(»« + *.)*■ (34) 

The mean anisotropy parameter A turns out to be 

A = ** + *? + <*? , (35) 

3Pk 2 (nlt + ki)n~ 2 
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The expansion 9 and shear scalar a 2 are given by 

a 3l 2 Qi 2 + Q2 2 + q 3 2 f ^ a , 

= — — , a = T . (36) 

rdt + ki 2k 2 (nlt + h)^ 

For stiff matter (w = 1), we have p = p. Thus the energy density and 
pressure of the universe become 

2np = 2ko = 5 7—. (37) 

{nlt + k 1 )™ +2 kinlt + k^- 

The entropy of universe is given by 
TdS = -[Qkl 3 (n - 2){nlt + h)^ 3 - U qi 2 + q 2 2 + q 3 2 )(nlt + h)-^ 1 ]. (38) 

K K 

Also, Eq.flHD leads to 

P = £, (39) 

where c is an integration constant. It is mentioned here that this value of p, 
when compared with the value obtained in Eq. (j3"Tj) . gives a constraint 

nc + 3kl 2 = (40) 

which holds only when n = 2 and q 2 + q 2 2 + q^ 2 = 0. The function of Ricci 
scalar, f(R) is 

f(R) = (nit + k x )^R + 3kl 2 (n - 2){nlt + h)^^ , (41) 

where R = R\ — 6l 2 (n — 2){nlt + ki)~ 2 . For a special case n = |, f(R) turns 
out to be 

/W =2 ( -9P ); -— <-9P» S (42) 

which gives /(-R) in terms of i?. 

3.4 Model of the Universe when n = 

Here the metric coefficients take the form 

A = Pl k 2 exp(lt)eM- qieX ^~ lt) ], (43) 

klk 2 

B = p 2 k 2 e^{lt)e X p[- q2e f l { - lt \ (44) 

C = p 3 k 2 exp(lt)ex P [- q3eX ^- U) }. (45) 
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The directional Hubble parameters Hi and the mean generalized Hubble 
parameter will become 

H i = l + ^-exp(-lt), H = l (46) 

kk 2 

while the volume scale factor turns out to be 

V = k 2 3 exp(3lt). (47) 

The mean anisotropy parameter A becomes 

A=[ ft a + ^ + fr a ]e,p(-2ft) (48) 

1 3Pk 2 k 2 2 J PV ; V ; 

while the quantities 6 and a 2 are given by 

9 = 31, a 2 = [^±^pl]eM-2lt). (49) 

For stiff matter, the energy density and pressure turn out to be 

2 Kp = 2 Kp = J"P^m _ (ft 3 + «■' + g^zM. (50) 

The corresponding entropy is 

TdS = -[-f2A;A; 2 / 3 exp(/t) - — (<?i 2 + <? 2 2 + g 3 2 ) exp(-Zt)]. (51) 

K /CK2 

The constraint equation with the condition, qi 2 + q 2 2 + q^ 2 = 0, is given by 

kc - 6kk 2 H 2 exp(4/t) = 0. (52) 

The function of Ricci scalar, f(R), takes the form 

f(R) = ^exp(-2lt){R - 12/ 2 ) (53) 

lK 2 

which reduces to 



U3/2 

f(R) = f-^[R-l2l 2 } (54) 

using the constraint Eg. (152)) . This corresponds to the general function f(R) 



f(R)=J2"nR n , (55) 

where n may take the values from negative or positive. 



4 Exact Bianchi Type V Solutions 

Here we shall find exact solutions of the Bianchi type V spacetime. 

4.1 Field Equations 

The metric for the Bianchi type V spacetime is 

ds 2 = dt 2 - A 2 (t)dx 2 - e 2mx [B 2 (t)dy 2 + C 2 (t)dz\ (56) 

where A, B and C are cosmic scale factors and m is an arbitrary constant. 
The corresponding Ricci scalar is 

A B C 3m 2 AB BC CA. 

R = - 2[ A + B + C-^ + AB + BC + CA^ (5?) 

With the help of Eq. fllj) , we can write 

B C 2m 2 AB CA AF F K , 
-B-C-^ + AB + CA + AF-F-F( p + p) = > (58) 

A C 2m 2 AB BC BF F k, 
-A-C-lP + AB + BC + BF-F-F^ + P >- - < M > 

A B 2m 2 BC CA CF F « , 
-A-~B-lV + BC + CA + CF-F-F^ + '' )=0 - (60) 

The 01-component can be written by using Eq.(j2J) in the following form 

ABC 
2 A-B-C=°- (61 > 

We discuss solution of these equations using the same procedure as for the 
Bianchi type I solutions. 

4.2 Solution of the Field Equations 

Here we get the same equations Eqs. fl21|) - fl23l as obtained previously with 
the difference of the constraint equations (using Eq.f loT]) ) 

pi = 1, p 2 = P3" 1 = P, Qi = 0, q 2 = -qz = Q- (62) 

Consequently, the metric functions become 

A = a, B = aPexp[QJ^}, C = aP' 1 exp[-Q J -^]. (63) 
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4.3 Model of the Universe when n/0 

The metric coefficients and the directional Hubble parameters are the same 
as given in pQ. Further, we note that the mean generalized Hubble parameter 
H, the volume scale factor V and expansion scalar 9 turn out to be the same 
as for the Bianchi type / spacetime while the value of shear scalar a 2 is 

a 2 = & (64) 

k 2 (nlt + ki)^ 
The mean anisotropy parameter for this model becomes 

A ^9. (65) 

3l 2 k 2 (nlt + h)^- 2 ' 

The energy density and pressure of the universe turn out to be 

3kl 2 Q 2 + k 2 m 2 ,_. 

np = no = = T . (bo) 

{nlt + h)^ +2 k{nlt + h)^ 
The entropy of the universe is given by 

1 9/ 

TdS = -[Qkl 3 (n - 2)(nlt + h)^ 3 - -(Q 2 + k 2 m 2 )(nlt + h)'" 1 }. (67) 
K k 

The constraint equation turns out to be same as for the Bianchi type I 
spacetime with the condition Q 2 + k 2 m 2 = 0. f(R) is given by 

f(R) = -(nit + h)^R + 3kl 2 (n - 2){nlt + k 1 ) = ^ 1 , (68) 

where R = 6l 2 (n — 2){nlt + ki)~ 2 + 8m 2 (nlt + ki)~^~]. For a special case, 
n — |, it follows that 

k-9l 2 ±V8ll i + 32m 2 R,„ 
f(R) = jl 2R ] R 



9kl 2 r -9/ 2 ± V81/ 4 + 32m? R, 



-3 



2 L 2R 

which gives f(R) in terms of R only. 
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4.4 Model of the Universe when n — 0. 

Here the metric coefficients, the mean generalized Hubble parameter, volume 
scale factor and expansion scalar are the same as given for Bianchi type / 
spacetime. However, the value of shear scalar a 2 is 

° 2 = W eXP( ~ 2/t) ' (69) 

and the mean anisotropy parameter A here becomes 

2q»exp(-2ft) 

3l 2 k 2 k 2 2 V ; 

The energy density and pressure of the universe are 



3fc/ 2 exp(-2/t) _ (Q 2 + m 2 k 2 ^ 
k 2 kk 2 



np = np = ~ - ,T J ' r J"'*"' ) exp(-4/t). (71) 



The entropy is 



1 2/ 

TdS = -\-12kk 2 l 3 exp(lt) - -—(Q 2 + m 2 k 2 )exp(-lt)]. (72) 

K kk 2 

The constraint equation with the condition, Q 2 + m 2 k 2 = 0, turns out to be 
same as given by Eq. fl52|) . The function of Ricci scalar, f(R), takes the form 

f(R) = _^ exp (_2/t)( J R - 12/ 2 ) (73) 



'2 
,2, 



2k 2 
with R= -12l 2 + (^r) exp(-2lt). Here we can get f(R) in terms of R 

/(i?)= 16^ [i?2 - 144/4] (74) 

which also corresponds to the general function f(R) given by Eq.(|55l). 



5 Summary and Conclusion 

This paper is to devoted to study the universe expansion in metric f(R) grav- 
ity. We have found exact solutions of the Bianchi types I and V spacetimes 
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using the non-vacuum field equations. These solutions correspond to two 
models of the universe, i.e., a singular model and a non-singular model. We 
have evaluated some important cosmological physical quantities for these 
solutions such as expansion scalar 9, shear scalar a 2 and mean anisotropy 
parameter A. The entropy of the universe is also found. We would like 
to mention here that solutions for both the spacetimes correspond to per- 
fect fluid [23 ES] in GR. It is found that the general function f(R) includes 
squared power of the Ricci scalar for the non-singular model. 

Firstly we investigate singular model (n ^ 0) of the universe which has 
singularity at t = t s = — ^y. The physical parameters Hi, H 2 , H%, H, 9, and 
a 2 are all infinite at this point for n > but the volume scale factor vanishes. 
The mean anisotropy parameter A is also infinite at this point for < n < 1 
and it will vanish for n > 1. The function of the Ricci scalar, energy density, 
pressure and T are also infinite while the metric functions vanish at this point 
of singularity. The model also suggests that the expansion and shear scalar 
decrease for n > with the passage of time. The mean anisotropy parameter 
also decreases for n > 1 with the increase in time. This indicates that after 
a large time the expansion will stop completely and the universe will achieve 

2 

isotropy. The isotropy condition, i.e., ^- — > as t — > 00, is also satisfied. 
The entropy of the universe is infinite for n > ~. Moreover, the second law 
of thermodynamics is valid for n > 2 with the conditions qi 2 + q 2 2 + q^ 2 = 
and Q 2 + k 2 m 2 = 0. Thus we can conclude from these observations that the 
model starts its expansion from zero volume with infinite energy density and 
pressure at t — t s and it continues to expand with time. 

The second model is non-singular model (n = 0) of the universe. The 
physical parameters Hi, Hi, H3, a 2 and A are all finite for all finite values 
of t. The mean generalized Hubble parameter H and expansion scalar 9 is 
constant while f(R) is also finite here. The metric functions do not vanish 
for this model. The entropy of the universe is finite while the second law of 
thermodynamics is valid for kk 2 < with the conditions qi 2 + q^ + q^ 2 = 
and Q 2 + k 2 m 2 = 0. The energy density and pressure become infinite as 
t — > — 00 which shows that the universe started its evolution in an infinite 
past with a strong pressure and energy density. The isotropy condition is 
also verified for this model. The volume scale factor increases exponentially 
with time which indicates that the universe starts its expansion with zero 
volume from infinite past. 

We would like to mention here that qualitative analysis of the solutions 
found has not been done. However, it would be interesting to perform simu- 
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lation procedure for the consistency of the results with Wilkinson Microwave 
Anisotropy Probe (WMAP) data as given in the literature. 

References 

[1] Sharif, M. and Shamir, M.F.: Class. Quantum Grav. 26(2009)235020. 

[2] Carmeli, M.: Commun. Theor. Phys. 5(1996)159; Riess, A.G. et al. 
(Supernova Search Team): Astron. J. 116(1998)1009; Bennett., CL. et 
al.: Astrophys. J. Suppl. 148(2003)1; Riess, A.G. et al.: Astrophys. J. 
607(2004)665; Spergel, D.N. et al.: Astrophys. J. Suppl. 170(2007)377. 

[3] Carroll, S.M.: Living Rev. Relativity 4(2001)1. 

[4] Nojiri, S. and Odintsov, S.D.: Int. J. Geom. Meth. Mod. Phys. 
115(2007)4. 

[5] Lobo, F.S.N.: The Dark Side of Gravity: Modified Theories of Gravity, 
invited chapter to appear in an edited collection " Dark Energy- Current 
Advances and Ideas"; larXiv: 0807. 1 640 . 

[6] Faraoni, V.: f(R) Gravity: Successes and Challenges, Presented 
at SIGRAV 2008, 18th Congress of the Italian Society of General 
Relativity and Gravitation, Cosenza, Italy September 22-25, 2008; 
larXiv:0810.2602l 

[7] Sotiriou, T.P and Faraoni, V.: Rev. Mod. Phys. 82(2010)451. 

[8] Nojiri, S. and Odintsov, S.D.: Problems of Modern Theoretical Physics, 
A Volume in honour of Prof. Buchbinder, I.L. On the occasion of his 
60th birthday, p.266-285, (TSPU publishing, Tomsk); ar Xiv:0807.0685l 

[9] Multamciki, T. and Vilja, I.: Phys. Rev. D74(2006)064022. 

[10] Multamaki, T. and Vilja, I.: Phys. Rev. D76(2007)064021. 

[11] Cognola, G., Elizalde, E., Nojiri, S., Odintsov, S.D. and Zerbini, S.: 
JCAP 0502(2005)010. 

[12] Capozziello, S., Stabile, A. and Troisi, A.: Class. Quantum Grav. 
24(2007)2153. 

14 



[13] Sharif, M. and Kausar, H.R.: Non-vacuum Static Spherically Symmetric 
Solutions in f(R) Gravity, submitted for publication. 

[14] Hollenstein, L. and Lobo, F.S.N.: Phys. Rev. D78(2008) 124007. 

[15] Azadi, A., Momeni, D. and Nouri-Zonoz, M.: Phys. Lett. 
B670(2008)210. 

[16] Momeni, D. and Gholizade, H.: Int. J. Mod. Phys. D18(2009)l. 

[17] Sharif, M. and Shamir, M.F.: Mod. Phys. Lett. A25(2010)1281. 

[18] Lorenz-Petzold, D.: Phys. Rev. D29(1984)2399. 

[19] Lorenz-Petzold, D.: Astrophys. Space Sci. 114(1985)277. 

[20] Berman, M.S.: Nuovo Cimento B74(1983)182. 

[21] Kumar, S. and Singh, C.P: Int. J. Theor. Phys. 47(2008)1722. 

[22] Singh, C.P: Brazilian J. Phys. 39(2009)4. 

[23] Paul, B.C., Debnath, PS. and Ghose, S.: Phys. Rev. D79(2009)083534. 

[24] Nojiri, S. and Odintsov, S.D.: Phys. Rev. D68(2003)123512. 

[25] Singh, C.P., Ram, S. and Zeyauddin, M.: Astrophys. Space Sci. 
315(2008)181. 

[26] Kumar, S. and Singh, C.P: Astrophys. Space Sci. 312(2007)57. 



15 



